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H ■ Abstract 

'^ ' We present analytic results for the current in a system moving in 

5 ' an arbitrary periodic potential and driven by weak Gaussian noise 

Qj , with an arbitrary power spectrum which are valid to order (tc/tr)'^, 

where tc is the largest characteristic time of the noise, and tr is the 

characteristic intrawell relaxation time. The dependence of the current 

on the shape of the potential, and on the shape of the power spectrum 

C^ ! of the noise is illustrated. It is demonstrated that the direction of the 

current is opposite when the power spectrum of the noise has minimum 
or maximum at zero frequency. A simple physical mechanism for this 
behavior is suggested. The behavior of the system in the limit of slow 
noise {tc » tr) is also discussed. 

Recently some intriguing work has appeared on the subject of stochas- 
tically driven ratchets, |jl], 0] and the transformation of noise in spatially- 
periodic (or phase-periodic) systems into a current. Such systems are noise 
rectifiers and thus are of great interest in relation to the understanding of 
machinery which operates in "the Brownian regime" [Q that is, on the very 
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small scale where fluctuations play a major role, and where the basic macro- 
scopic methods of controlling energy flow no longer remain valid. The most 
immediate examples are biological systems, [^ but these ideas may not be 
irrelevant in the area of applied technology, where a great interest has devel- 
oped recently in the possible construction of nanoscale devices. P] 

A simple model for a ratchet is a noise-driven overdamped nonlinear dy- 
namical system described by the stochastic differential equation 

i = -U\x) + f{t), U{x) = U{x + X), (/(t)) = 0, (1) 

where U{x) is a periodic potential, such as the one illustrated in Fig. 1, and 
/(t) is zero-mean noise of some type. The theoretical problem is to find the 
stationary current density j = {x{t)) in the ratchet given the shape of U{x) 
and the properties of the noise f{t), and then to find the most appropriate 
conditions for the transformation of the noise into the current. If the noise is 
white then the system is (quasi)thermal, the second law of thermodynamics 
applies, and j = 0. If the noise is not white, i.e., for colored noise, the system 
is no longer in thermal equilibrium, and in general j 7^ 0. Since onset of a 
current means breaking the "right-left" symmetry, currents may only arise, 
in the case of additive noise, if the potential U{x) is asymmetric with respect 
to its extrema. The onset of a current can be viewed also as an example 
of "temporal order coming out of disorder" , since the current is apparently 
time-irreversible, whereas stationary noise does not distinguish "future" from 
the "past"; we notice, however, that Eq. (1) implies relaxation and is thus 
time-irreversible itself. 

We consider the case where f{t) is zero-mean Gaussian noise with a 
frequency-dependent power spectrum 

/oo 
dt expM)0(t), 0(t) =< /(t)/(0) > . (2) 
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If the characteristic noise intensity D = max$(u;) is small then most of the 
time the system performs small- amplitude fluctuations about the minima of 
the potential. Occasionally it "jumps" from the minimum it occupied to the 
one on the right or left, with the probabilities per unit time W+ and W^, 
respectively. These jumps give rise to the current j = A(W/'+ — W-). The 
dependence of the transition probabilities W± on the noise intensity D is of 
the activation type for Gaussian noise, W± = const, x exp{—R±/D). The 



characteristic activation energies R± of the escapes from a given minimum 
Xo of the potential over the right (x+) or left (x_) maximum are given by 
the solution of a variational problem formulated in |Q. In the case where the 
bandwidth of the spectrum $(c(j) greatly exceeds the reciprocal relaxation 
time of the system t~^ = U"{xq) the activation energies are R± = 2F{0)AU+ 
7± -^"(0), where 



AU = U{b) - U{a) > 0, 7± = T* dx U'{x)[U"{x)f > 0, (3) 
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FM=D/$H, n^) = ^^> |F"(0)/F(0)|«t^ (4) 

The current is then given by 

J = Al^^[exp(-7+F"(0)/D) - exp(-7_F"(0)/D)], (5) 

where Wk is the Kramer's activation rate 

Wk = ^D/^i^) ^ exp(-2F(0)AU/D). (6) 

^U"{xo)\U"{x^)\ 

Eq. 5 is a principle result of this letter. It is immediately obvious from (5) 
that: (i) the noise color does give rise to the onset of a current due to the fluc- 
tuational interwell transitions, and (ii) the direction of the current depends 
crucially on the shape of the spectral density since F"{0) can take both posi- 
tive or negative signs. We emphasize that, although the corrections 7±-F"(0) 
to R± are small compared to the main term, they are not small compared to 
the noise intensity D and can change W± by orders of magnitude (we have 
neglected the corrections to the prefactor in W± due to the noise color and 
used the standard Kramers expression for this prefactor valid for white-noise 
driven systems). In fact, except the special case where U{x) is symmetric 
with respect to Xq, the ratio of the probabilities W^/W^ is exponentially large 
or small for small intensity of the colored noise, and therefore the transitions 
in one direction dominate overwhelmingly over the transitions in the opposite 
direction, so that \j\ ^ XWy where Wy = ina.x{W+,W-). The direction of 
the current is determined by the interplay of the shape of the potential and 
the features (the shape of the power spectrum, in the present case) of the 
noise. 



The dependence of the current on the noise color is different depending on 
the sign of F"{0). When F"{0) > there is a saturation effect where, given 
the noise strength D and the shape of the potential, the current is maximal 
for F"{0) given by 

(7- - 7+) 

Thus, for F"(0) > 0, F^{0) is the optimal noise color for a given noise 
strength and ratchet potential. 

Eq. (5) also provides an answer to the following question: given the barrier 
height AU, what is the shape of the well for which the current (5) will be most 
pronounced for weak noise color, that is, what is the most effective shape of 
the ratchet? This shape is given by the extreme value of 7± for given AU and 
period A. It is straightforward to show that the corresponding variational 
problem does not have a differentiable solution: the single-valued potential 
has to be of the shape approaching the one pictured in Fig. 2. A simple 
analysis can be done, e.g., for U'{x) = Ai + 742tanh[(a; — xo)/a], —A/2 < 
X < A/2 {Ai + A2tanh[(xo — xo)/a] = 0): for small a/A the shape of the 
corresponding potential is close to that of the sawtooth. 

The above results apply immediately to the most often studied case where 
f(t) is exponentially correlated, 

0(t) = -^exp(-|t|/r), $H =D/(1 + .;V). (8) 

The problem of fluctuational transitions induced by this noise has been in- 
vestigated for small noise intensities in very much detail (see 0, 0, 
and also |TD[ for a review). The noise (8) has one correlation time, r, and 



F"[0) = 2r^ > 0. Obviously, j vanishes to first order in r, in agreement with 
the results of [2]. It has a maximum as a function of r for r given by (7), 
and falls down for large 7±r^/D. The saturation of j vs. r is illustrated in 
Fig. 3. 

Naturally occurring noise will generally not be exponentially correlated, 
a situation which any realistic physical theory must accommodate. This is 
clearly the regime in which j ^ 0. The advantage of the result presented 
here is that it not only gives an analytic result to second order in tc/tr, but 
is valid for any Gaussian noise. A more general situation than (8) can be 



modeled ||T0[ by the noise with the power spectrum 






[uj^ - ul)^ + AV^u^'- 



where F is a measure of the bandwidth, and c^o the frequency of the noise. 
This can be thought of as the power spectrum of the system, 

<mi{s)>=AVb5{t~s). (10) 

In this case F(0) = ujID/AVID, and F"(0) = (21^ - ujl)D/VJD. Obviously, 
F"(0) > when uJq < 2T^ and thus the maximum of the power spectrum (8) 
is at tu = 0. In this case the direction of the current for F ^ t~^ is the same 
as for exponentially correlated noise and the analysis given above directly 
applies. 

A completely different situation occurs if uJq > 2F^ and the power spec- 
trum (8) has a minimum at cj = 0. In this case F"{0) < 0. It follows 
from (5) that the direction of the current is opposite to that arising for ex- 
ponentially correlated noise. In contrast to the case F"{0) > considered 
above, in the present case increase in |F"(0)| does not give rise to saturation 
and then to decrease in j: the current is increasing exponentially with the 
increasing |F"(0)| where the approximation (3) is applicable. The current 
reversals have been found recently in numerical experiments, and in certain 
specific exactly solvable cases, by Doering and Horsthemke^. No plausible 
physical mechanism has been suggested as an explanation in |0. Not only 
does the present analysis give an analytic criterion for the current reversal, 
$"(0) oc — F"(0) > 0, but it also suggests a direct physical interpretation. 

The correction to the activation energy of a transition R in the case of 
"weakly colored" noise is due to the fact that it is not only the total work 
the noise does on the system on the way from the potential minimum to the 
barrier top (along the optimal path of the escape) that counts, as in the case 
of white noise, but the characteristic strength and duration of the pulse that 
gives rise to the escape. The parameters 7± just characterize the ratio of the 
squared value of the force to the duration of the pulse (cf. 0). The shorter 
the pulse the higher are the frequencies of the noise components involved. 
If the power spectrum of the noise $(a;) decreases with a; {F"{0) > 0) then 



the higher-frequency noise components are weaker on the average, and the 
escape probabihty decreases with the decreasing duration of the pulse. For 
$(u;) increasing with uj (-F"(0) < 0) the result is exactly opposite. The 
characteristic duration of the pulse is determined by the ratio of the distance 
between the potential minimum and the barrier top to the characteristic 
velocity [/', and for a given height of the potential barrier it scales as ([/')~^. 

Suppose we have a potential more steep to the left from the minimum, as 
shown in Fig.l. It is clear from (3) that in this case 7_ > 7+. The pulse which 
gives rise to the escape over the left barrier is shorter, and the characteristic 
frequencies involved are higher. Therefore, for F"(0) > (as in the case of 
exponentially correlated noise), the escape over the left barrier is less likely 
to happen than that over the right one, and the current flows to the right, 
whereas for F"(0) < it flows to the left. 

The behavior of the system when driven by slow noise itc » tr) is also 
of interest. The activation energies in this limit can also be calculated via 
the method of Ref. p, and are given by R± = D[U^ _|_]^/2(^(0) where It/^ ±1 
are the maximum values of \U'\ on the intervals {xo,x±). We note here in 
passing that there are nonanalytic [{tr/tcY^^] corrections to R± @, Q. The 
current which arises in the case of slow noise is then 

J = A[C+exp(-[f/;,jV20(O)) -C_exp(-[tZ;„,_]V20(O))]. (11) 

(C-t are the constants that allow for the prefactors in the expressions for 
W±; obviously, 0(0) is just the mean-square value of the noise). What is 
happening here also has a clear physical interpretation (cf. P]). Since the 
noise has an extremely long correlation time, in the range where U"{x) > 
the particle simply follows the force adiabatically according to U'{x(t)) = 
f{t). The fluctuation large enough to allow the particle to escape over a 
barrier is just the one that overcomes the restoring force —U'{x) for all x, 
and the probability of such a fluctuation is just W oc exp{—R/D). The 
current in this case is always in the positive direction for [f/^_]^ > [f/^_,_]^. 
There is again a saturation effect, and the current is maximized for 0~^(O) ~ 

4Din(ic„f/;,_/c+f/;„,j)/([f/;,_]2 - [u:^^^n 

We emphasize that the onset of current in the system considered is the 
result of it being away from thermal equilibrium. This means that dissipation 
and fluctuations are not interrelated via fluctuation-dissipation theorem. In 
the particular case considered dissipation was not retarded (the friction force 



is determined by the instantaneous value of the coordinate), and therefore 
the effect arises when the power spectrum of the noise displays dispersion (for 
Gaussian noise). A current can also arise in the situation of a white-noise 
driven system where dissipation is retarded. 

Lastly we remark that the system can be expected to exhibit behavior 
similar to stochastic resonance as the noise strength D is varied, with j 
exhibiting a maximum for some noise strength. Obviously for very small D 
the current is an increasing function of D. As the noise strength is increased 
the fluctuations will begin to wash out the effects due to the shape of the 
potential, leading to a decrease of the current for large D. 
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Figure 1: Typical ratchet potential U{x) 



Figure 2: The sawtooth potential which is the "optimal" limit of U{x) which 
gives rise to the greatest current. 
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Figure 3: Current vs. r for varying potential shape showing that the current is 
maximized for a specific r. 



